Azimuthal asymmetries for hadron distributions inside a jet in hadronic collisions 
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Using a generalized parton model approach including spin and intrinsic parton motion effects, and 
assuming the validity of factorization for large-pr jet production in hadronic collisions, we study the 
azimuthal distribution around the jet axis of leading unpolarized or (pseudo)scalar hadrons, namely 
pions, produced in the jet fragmentation process. We identify the observable leading-twist azimuthal 
asymmetries for the unpolarized and single-polarized case related to quark and gluon-originated jets. 
We account for all physically allowed combinations of the transverse momentum dependent (TMD) 
parton distribution and fragmentation functions, with special attention to the Sivers, Boer-Mulders, 
and transversity quark distributions, and to the Collins fragmentation function for quarks (and to 
! >— , the analogous functions for gluons). For each of these effects we evaluate, at central and forward 

rapidities and for kinematical configurations accessible at BNL-RHIC, the corresponding potentially 
maximized asymmetry (for tv + production), obtained by saturating natural positivity bounds (and 
the Soffer bound for transversity) for the distribution and fragmentation functions involved and 
summing additively all partonic contributions. We then estimate, for both neutral and charged 
pions, the asymmetries involving TMD functions for which parameterizations are available. We also 
study the role of the different mechanisms, and the corresponding transverse single spin asymmetries, 
for large-pT inclusive jet production. 

PACS numbers: 13.88.+e, 12.38.Bx, 13.85.Ni, 13.87.Fh 

in: 

I. INTRODUCTION 

Transverse single-spin and azimuthal asymmetries in high-energy hadronic reactions have raised a lot of interest in 
the last years (see e.g. Refs. [HQ and references therein). Huge spin asymmetries have been measured in the inclusive 
forward production of pions in high-energy pp collisions at moderately large transverse momentum. The general 
trend of the early pioneer measurements of the E704 Collaboration at Fermilab [U, 0| has been recently confirmed at 
much larger center of mass (cm.) energies at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
• ■ Laboratories (BNL) in similar kinematical configurations @. A surprisingly large transverse polarization of A 
,— ' \ hyperons produced in the forward region was also measured in unpolarized pp, pN fixed-target experiments, see e.g. 
Ref. 0- Also in this case it will be hopefully possible in the near future to check if this intriguing effect survives 
at the much larger energies reachable at RHIC and at the Large Hadron Collider (LHC) at CERN. Similar effects, 
• • . leading to azimuthal asymmetries both in the polarized and unpolarized case, have been measured in Drell-Yan (DY) 
processes in semi-inclusive deeply inelastic scattering (SIDIS) [If3 - ll3l |. and in hadron-pair production in e + e~ 

collisions |14l. Ha|. 

These results cannot be explained at leading-twist (LT) approximation in the usual collinear approach of pertur- 
bative QCD (pQCD), based on factorization theorems, to inclusive particle production in hadronic collisions. Here 
collinear means that intrinsic parton motion is neglected in the hard scattering processes and integrated over up 
to the large energy scale in the soft functions involved. On the contrary, at least in the kinematical regimes under 
consideration at RHIC, collinear next-to-leading order (NLO) pQCD gives a fair account of unpolarized cross sections 
(see e.g. Refs. [HE3]). 

Two different main theoretical approaches have been proposed in the framework of perturbative QCD in order to 
account for these measurements. One is the so-called twist-three collinear approach, which generalizes the leading 
order (LO) collinear framework with the inclusion of higher-twist quark-gluon correlations [18l420l |. This involves a 
new class of universal nonperturbative twist-three quark-gluon distribution and fragmentation functions that need to 
be modeled by fitting experimental data. Another formalism, which will be adopted in this paper, is the so-called 
transverse momentum dependent (TMD) approach, which takes into account spin and intrinsic parton motion effects. 
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Although the large single spin asymmetries (SSAs) of interest here were originally observed in single inclusive particle 
production in hadronic collisions, it is now clear that from the theoretical point of view these are not the cleanest 
processes to consider. First of all, these SSAs are twist-three effects in a series expansion in inverse powers of the large 
energy scale (here, the transverse momentum of the observed single hadron or jet). Several competing mechanisms 
can therefore play a role and mix up. Moreover, in the TMD formalism factorization has not yet been proved for these 
processes and its validity is much under debate presently, see e.g. Refs. [2ll - l26l | for inclusive two-hadron production. 
Factorization breaking would in turn imply non universality of the soft TMD distribution and fragmentation functions 
required to explain data. In the case of single particle production in hadronic collisions, therefore, the TMD approach 
can be seen at present as a useful generalization of the parton model where factorization is assumed as a reasonable 
starting point to be carefully scrutinized by comparison with available experimental results. 

As mentioned above, similar spin and azimuthal asymmetries were subsequently observed in SIDIS and DY pro- 
cesses, where two energy scales are involved, the required large energy scale (the momentum transfer Q in SIDIS, the 
lepton-pair invariant mass, M , in DY) allowing use of pQCD, and a small energy scale sensitive to intrinsic parton 
motion (the transverse momentum respectively of the final lepton pair in DY and of the produced hadron in SIDIS). 
For these classes of processes both the twist-three collinear formalism and a full gauge invariant TMD approach have 
been developed and factorization has been proven Moreover, it has been shown that when the value of the 

small observed scale is intermediate between the typical QCD nonperturbative scale and the large factorization scale, 
the two approaches are mutually consistent [3lM33| . 

In the TMD approach to DY(SIDIS) processes color gauge invariance is ensured by the inclusion of gauge links 
(Wilson lines), accounting for initial(final) state interactions among the struck partons involved in the hard process and 
the remnants of the parent hadrons (additional final state interactions are also present in the fragmentation process). 
Single-spin and azimuthal asymmetries are generated by TMD polarized partonic distribution and fragmentation 
functions, among which the most relevant from a phenomenological point of view are the Sivers distribution [34], 
and, for transversely polarized quarks, the Boer-Mulders distribution [36j and the Collins fragmentation function 
(similar functions can be defined for linearly polarized gluons, see e.g. Ref. [Hj]). 

For inclusive forward pion production the large transverse single spin asymmetry observed can be generated both 
by the Sivers and the Collins effects; unfortunately these contributions cannot be disentangled and one has to consider 
alternative measurements in order to separate the different mechanisms. This is at variance with the case of SIDIS 
and DY processes, where the Sivers and Collins effects (and several other possible contributions to the azimuthal 
asymmetries) can be disentangled. In hadronic collisions one has to resort to different processes, like e.g. the DY 
process (no fragmentation), single photon or jet production, two particle(jet) production with transverse momentum 
imbalance and so on. 

From this point of view, a very interesting process is pp — > jct + 7r+Y, where one observes a large pt jet and looks for 
the azimuthal distribution of leading pions inside the jet. In this case, one should observe a symmetric pion distribution 
for the fragmentation of an unpolarized parton jet, and a cos<f> (cos 20), distribution for a transversely (linearly) 
polarized quark(gluon) parton jet (</> indicates the azimuthal angle of the leading pion distribution around the jet 
axis). Therefore, despite the complexity of the measurement (which is however at reach and presently under active 
investigation at RHIC), this process might offer plenty of new information as compared to single inclusive pion 
production. It would in principle allow us to disentangle the contributions coming from the Sivers and the Collins 
effects. Other contributions involving different combinations of TMD distribution and fragmentation functions could 
also be disentangled. Finally, it could also help in identifying jets coming from quark or gluon fragmentation, since 
the pion azimuthal distribution is different in the two cases. At RHIC kinematics a careful tuning of the kinematical 
configuration considered can help from this point of view. 

Motivated by these considerations, in this paper we will present, in the approach of the TMD generalized parton 
model, and allowing for intrinsic parton motion both in the initial colliding hadrons and in the fragmentation process 
(which is crucial), the general expression for the polarized cross section for the process p^p — > jet + it + X, and 
the structure of the azimuthal asymmetries that can be measured in the distribution of leading pions around the 
jet thrust axis (coinciding in our scheme with the final scattered parton direction of motion). A very preliminary 
version of this study was first presented in Ref. 39] . A similar analysis was discussed in Ref. [4(| , which however 
considered intrinsic parton motion only in the fragmentation process, drastically reducing the possible contributions 
to the asymmetry. Indeed, in that case, only the Collins effect for quarks is at work. In fact, Ref. [4(| aimed at 
studying only the Collins fragmentation function (FF), which should be universal, in a more simplified theoretical 
scheme for which factorization has been proven. Our approach is different in some respects. It is more general and has 
in principle a richer structure in the observable azimuthal asymmetries, since intrinsic motion is taken into account 
in the initial hadrons also. However, since factorization has not been proven in this case, but is rather taken as a 
reasonable phenomenological assumption, the validity of the scheme and the universality of the TMD distributions 
involved require an even more severe scrutiny by comparison with experimental results. On the other hand, at the 
present theoretical and experimental stage, we believe that combined phenomenological tests of different approaches 



3 



are required to clarify the validity of factorization and, related to this, the relevance of possible universality-breaking 
terms for the TMD distributions. 

The plan of the paper is the following. In Sect. |n] we will summarize the TMD generalized parton model approach, 
which has been presented and discussed at length in a series of papers, see e.g. Refs. [38, 4l|, [42j. We will then present 
the expression of the polarized cross section for the process of interest, discussing in detail the different partonic 
contributions to the process; we will finally list the azimuthal asymmetries that can be measured and their physical 
content. In Sect. UTTl we will present phenomenological results for the azimuthal asymmetries discussed in the kine- 
matical configuration of the RHIC experiments, at different cm. energies and for central and forward rapidity jet 
production. In particular, we will first present results for the totally maximized effects, by taking all TMD functions 
saturated to natural positivity bounds and adding in sign all possible partonic contributions. This will assess the 
potential phenomenological relevance of each effect. We will then consider more carefully those effects involving the 
Sivers and Boer-Mulders distributions and the Collins fragmentation function, for which phenomenological parame- 
terizations obtained by fitting combined data for azimuthal asymmetries in SIDIS, Drell-Yan and e + e~ collisions are 
available. Section llVl contains our final remarks and conclusions. 



II. FORMALISM 



In this section we present and summarize the expressions of the polarized cross section and of the measurable 
azimuthal asymmetries for the process A^B — » jet + tt + X, where A and B are typically a pp or pp pair. Since most 
of the formalism has been already presented in Refs. [HI, |4l|, [42[, we will shortly recall the main ingredients of the 
approach, discussing more extensively only relevant details specific to the process considered. 

Within a generalized TMD parton model approach including spin and intrinsic parton motion effects, and assuming 
factorization, the invariant differential cross section for the process A(Sa) B — > jet + tx + X can be written, at leading 
twist in the soft TMD functions, as follows: 

E da A(S A ) B^ c t+«+X f dXadxb 

d 3 Pj dzd^ = £ J 16A aXb s dk ^ d k ^<K f«/A,sAx a ,^ a ) P < bK f b/B ( Xb ,k^) 

a,b,c,d,{\} 

X ^A e ,A di A a) A 6 ^A il X (|i A;,Ai*(« + * + «) 6 5 1! ,A i (« > fc±>r)- W 

In a LO pQCD approach the scattered parton c in the hard elementary process ab — > cd is identified with the 
observed fragmentation jet. Let us summarize briefly the physical meaning of the terms in Eq. ([1]). Full details and 
technical aspects can be found in Refs. |38L EH. l42j| . 

We sum over all allowed partonic processes contributing to the physical process observed. {A} stays for a sum over 
all partonic helicities, A = ±1/2 (±1) for quark(gluon) partons respectively. x a ,b and fcj_ a ,b are respectively the initial 
parton light-cone momentum fractions and intrinsic transverse momenta. Analogously, z and kj_^ are the light-cone 
momentum fraction and the transverse momentum of the observed pion inside the jet with respect to (w.r.t.) the jet 
(parton c) direction of motion. 

P\^\' SA fa/A,SA i x ai kj_a) contains all information on the polarization state of the initial parton a, which depends 
in turn on the (experimentally fixed) parent hadron A polarization state and on the soft, nonperturbative dynamics 
encoded in the eight leading-twist polarized and transverse momentum dependent parton distribution functions, which 
will be discussed in the following. P\^\/ Sa is the helicity density matrix of parton a. Analogously, the polarization 

state of parton b inside the unpolarized hadron B is encoded into p^^, fb/B( x b, &_Lb)- 

The ^A c ,A.;A a ,A 's are the pQCD leading-order helicity scattering amplitudes for the hard partonic process ab — s- cd. 

The y (z, k± n )'s are the soft leading-twist TMD fragmentation functions describing the fragmentation process 
of the scattered (polarized) parton c into the final leading pion inside the jet. 

As already said, we will consider as initial particles A, B, two spin 1/2 hadrons (typically, two protons) with hadron 
B unpolarized and hadron A in a pure transverse spin state denoted by Sa, with polarization (pseudo) vector P A . 

Ej and Pj are respectively the energy and three-momentum of the observed jet. 

Unless otherwise stated, we will always work in the AB hadronic cm. frame, with hadron A moving along the 
+Z cm direction; we will define (XZ) cm as the production plane containing the colliding beams and the observed jet, 
with (pj)x cm > 0. We therefore have, neglecting all masses (see also Fig. [T]): 
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FIG. 1: (color online) . Kinematical configuration for the process A(Sa)B — > jet + tv + X in the hadronic cm. reference frame. 
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J>c = pf = ( E },P}T,0,pj L ) = Ej(l, sin 0j,O,cos0j) = t (cosh ?7j , 1, 0, sinhjfi) 

= i?,r ( 1, sin # w cos </>,r, sin # T sin </>,,•, cos 6*^ ) , (2) 

where jj is the jet (pseudo)rapidity, ryj = — log[tan(#j/2)]. 

Notice that, since the observed jet is identified with the scattered parton c, the helicity frame of the fragmenting 
parton, whose z axis, Zj, is along the direction of motion of parton c, is related to the hadronic cm. frame by a simple 
rotation by 9j around Y cm = y y In this frame (zj identifies also the jet light-cone direction) we have: 

= # = -$(1,0,0,1) 



Pi = E„ k ±n , J El - k\ \ = E„, k ±7r cos^f, fc x * sin^f, J E% - k 2 Xn , (3) 



where 0^ is the azimuthal angle of the pion three-momentum around the jet direction of motion, as measured in the 
fragmenting parton helicity frame. 

The light-cone momentum fraction of the observed pion is given by 



_ p+ = p+ _ E n + y/E%-Kj„ 

pt'pf 2E } ■ W 
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We can also write, respectively in the jet (parton c) helicity frame and in the hadronic cm. frame: 



Pn = k ±n cos <j>% xj + k ±7r sin <j>% y. + J El - k\ z 



j 



fc ±7r cos <^>f cos 6j + J El- fc^ sin 6 3 X cm + k ±7T sin (f)"Y cm (5) 



- fc ±7r cos <f>" sin 6»j + sjE*- k 2 j_ w cos 0j Z cm . 

Let us stress that in our notation intrinsic transverse momenta, i — a, b, tt, are always three- vectors and 

fc_i_j = |fej_i|. This has to be kept in mind when comparing with literature, where often intrinsic momenta are intended 
as 4- vectors and fc^ = — fe^. From Eq. ((SJ) it is easy to see that the pion intrinsic transverse momentum is given, in 
the hadronic cm. frame, by 

k±TT = fc ±7r cos0f cos6»jX cm + k ±7T sin (f>"Y cm - fc i7r cos</)f sin6>jZ cm . (6) 

Defining by 4>k the azimuthal angle of the pion intrinsic transverse momentum, k^, as measured in the hadronic 
cm. frame, see Fig. [TJ from Eq. ([6]) we easily see that 

tan (A? . , 

tan<^ = 7 

COS#j 

Notice that, for central-rapidity jets (0j = tt/2), <pk — n/2. Therefore, azimuthal asymmetries modulated in terms 
of (pk are artificially suppressed in the central rapidity region, while the physically relevant angle is <fi^ . Instead, in 
the forward rapidity region, when cos#j — > 1, the two angles are practically coincident. Notice also that in Ref. p0| . 
where only the forward rapidity region was considered, the angle <pk (called <ph there) was adopted. 

Let us now come back to the soft TMD partonic distribution and fragmentation functions entering the differential 
cross section for the process A(Sa) B — > jet + tt + X, Eq. ([T]). Consider first the polarized soft process at the 
distribution level, A(Sa) ~ > a + X; as said, in Eq. (|T|) all information on this process is encoded in the factor 
pO/A,SAj!^ g A (x a , k± a ). This factor depends on the polarization state (fixed by experimental conditions) of the 



parent hadron A, described by its own helicity density matrix p A ' SA , and on generalized soft distribution functions 



for the process A(Sa) — > a + X, F\"\? (x a , k± a ) 



The functions F^ a '^J are related to the well-known leading- twist hand-bag diagram for deeply inelastic scattering. 

Analogous relations hold for parton b inside the (un)polarized hadron B. 

Rotational invariance and parity conservation for strong interactions imply some very general relations for these 
nonperturbative functions: 

F^(x a ,k ±a ) = J£&(*.,*x.)e^- A *>*- , (9) 

where <fi a is the azimuthal angle of parton a intrinsic transverse momentum, k± a , in the parent hadron A helicity 
frame (coinciding in the present case with the hadronic cm. frame). Notice that the reduced soft functions F on the 
right-hand side (without the hat) do not depend anymore on azimuthal phases. Moreover, 

F:^i(x a ,k ±a ) = (-i)^---)(-i)( A A-A )+(Ak-vj F A.x (a . o)A . ±a)i (10) 

where Sa and s a are the spins of the parent hadron and of the parton respectively. Notice that for spin 1/2 colliding 
hadrons, Sa — 1/2, the factor (— l)^ - ^) is positive(negative) for quark(gluon) partons, and therefore some parity 
properties of the soft functions are different for quarks and gluons. By combining complex conjugation and parity 
properties, one can show that some of these functions are purely real or purely imaginary. As a result, for spin 1/2 
hadrons only eight independent soft distributions survive at leading twist. They can be easily related to the TMD 
distributions widely discussed in the literature. Before listing them for completeness, let us recall another important 

A A' 

property of the hadronic functions F x a .? [x a , k± a ), coming from total angular momentum conservation in the forward 
direction, that is for k± a — > 0: 

A A' /fcxaV AA_Aa ~ (A ' A ~ Vft)l ~A A' 

F A",'v"( a! »» fc -J-«»)~ ( -Jf) ^.A^a.fcxa), (11) 
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where M is a typical hadronic mass scale and the F's stay for the remaining part of the functions that, depending on 
the details of dynamics, may or may not vanish in the collinear configuration. The relevant hadronic functions and 
their connection with the leading-twist TMD distributions is therefore, in the quark case (for clarity we adopt here 
both the notation of Ref. [38[ and that of the Amsterdam group [2l|, HH, HH) 



-+ 
.++ 



F+- 
F~ + 



ReF q ++ 



ReF q : 



im^; 



X, kj_ 

x,k ± 
x, k± 

x, k± 

x, k±_ 

x, k±_ 

x, k±_ 

x, k 



f q /A(x,kj_) = ff(x,k ± ) 
^hf q /A{x, k±) = g\ L {x, k ± ) 
h\(x, k±) 
k\ j_ 

±A N f q/A ,(x,k ± ) = 
^h±Z(x,k ± ) 
-lA N f qt/A (x,k ± ) 



k 

Ym 



7T77 fir ( x ' k±) 



—— h-, q (x, k±) 
2M 1 v ' ' 



(12) 



Analogous relations hold for gluon partons with the changes discussed above, due to the different spin of the parton 
and leading to different parity properties. Notice that instead of transversely polarized quarks we will have linearly 
polarized gluons. Of course, the same relations hold also for the B — > b + X process. However, this time the Xb and 
Zg axes of hadron B helicity frame are opposite to those of the hadronic cm. frame. 

Concerning the fragmentation process, since here we are considering only pions (in general, unpolarized hadrons), 
the discussion of the soft fragmentation functions, D x x , , is much simplified. In practice, only two independent TMD 
fragmentation functions survive: one with diagonal parton helicity indexes, related to the TMD unpolarized FF, 



D±±{z, k ±n ) = -D±±(^, fc_i_7r) = D w/c (z, k ±7r ) , 



(13) 



and a second one with off-diagonal parton helicity indexes, D+_ (z, k±^). This second function is purely imaginary 
for quark partons and is related to the well-known Collins function [37[ , describing the fragmentation of a transversely 
polarized quark into a non collinear unpolarized hadron. Instead, for gluon partons the analogous Collins-like function 
is purely real and is related to the fragmentation of linearly polarized gluons again into an unpolarized hadron. 

It is very important to realize that these off-diagonal quark and gluon TMD FFs have different behaviours as a 
function of the azimuthal angle of the observed pion around the direction of motion of the fragmentation jet. More 
specifically, for quarks we have: 



&£(z,k x «)=±Dl'*(z,kx« 



■,±i4 



while for gluons the analogous relation reads 



D n J4(z, k ±n ) = D^{z, k^)e ±a ^ 



(14) 



(15) 



As mentioned above, the azimuthal independent parts of these TMD FFs are related to the probability for a trans- 
versely (linearly) polarized quark(gluon) of fragmenting into a non collinear unpolarized hadron, the Collins (Collins-like 
for gluons) fragmentation function (38j , 



Dl^(z,k ± „) = jA% t ( z ,M = ^ ff i( z -^) 



D^(z.k^) = ^A»D„ /77 {z,k ± J 



N; 



(16) 



Again, total angular momentum conservation in the forward direction dictates the power behaviour of the TMD 
FFs for k±Tr -> 0: 



Dl x ,(z,k^) 



M 



A c~ A cl 



Dl x ,(z,k ±v ). 



(17) 
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Also in this case, the behaviour of the off-diagonal FFs is different for quarks and gluons. 

Let us finally comment on the helicity amplitudes for the partonic hard scattering processes entering Eq. ([T]). Again, 
details have been already presented in Refs. [H, [42| and we limit here to summarize some useful properties. Due to 
intrinsic partonic motion in the distributions and in the fragmentation process, the general kinematical configuration 
for the partonic process ab — > cd is not planar in the hadronic cm. frame. Since partons can be in general polarized in 
the process, azimuthal phases are therefore essential and must be properly taken into account. For massless partons, 
due to helicity conservation in the quark-gluon vertex and parity invariance only three independent helicity amplitudes 
survive: 





-+ = M*__._ 


= M°e itfn 






+ = H-A 


= M°e lV2 


(18) 




- = 


_ + = M 3 °e^ 3 . 





Here ± stays for A = ±1/2 for quarks and A = ±1 for gluons. Mf {i = 1,2,3) are the three independent helicity 
amplitudes in the canonical partonic cm. frame, that is a frame where partons a, b move along the ±z direction 
respectively and the scattering plane coincides with the (xz) plane. The phases ifi collect all azimuthal phases coming 
from rotations and boosts connecting the canonical partonic cm. frame with the hadronic cm. frame adopted in the 
paper. Their general expression is rather involved. All details and the explicit expressions of the Mf and <pi can 
be found in Refs. [H, H2|. Here we limit ourselves to notice that for parton c lying in the (XZ) cm plane, as in the 
present case, the phases ipi are odd under fe^ a .b — fc_La.b- This property, as we will see below, is very helpful in 
selecting physically observable effects out of the many contributions present in Eq. ([1]) because of the nonplanarity 
of the partonic process, which however do not survive at the hadronic level under integration over intrinsic parton 
momenta. 

We now concentrate on the partonic kernels entering the expression of the polarized cross section, Eq. (|TJ) : 

{A} 

X MA e ,A <l iA al A i M^, A)li A;,Ai^A e) A4(*.*'X*). (19) 

One has to evaluate the kernels for each of the eight distinct partonic channels contributing to the cross section, 

qq -> qq, qg -> qg, qg -> gq, gq -> qg, gq -> gq (20) 
99 -> 99, 99 -> 99, 99 ~> 99 , 

where in the first line q stays for both quarks and antiquarks in all allowed combinations. 

In practice, the calculation is performed by summing explicitly over all helicity indexes and inserting the appropriate 
expressions for the helicity density matrices of partons a, b and for the polarized distribution and fragmentation 
functions, as detailed above. Furthermore, after factorizing explicitly all azimuthal dependences, including those 
coming from the hard-scattering helicity amplitudes, collecting them and using symmetry properties under k± a ,b — > 
—k± a ,b, one gets the final expression for the kernels, containing only physically allowed terms at the hadronic level. 

We will not present explicitly the kernels for all channels. Instead, we limit ourselves to give the kernels for 
the qq — > qq and gg — > gg channels, which contain the maximal number of terms and give examples of possible 
contributions involving both quark and gluon distribution and fragmentation functions. Moreover, we will directly 
present the combination of kernels, E(0s 4 ) ± £(</>s A + n) entering the numerator and the denominator of the single 
spin azimuthal asymmetries discussed in the sequel. We will also omit for shortness the explicit dependences on the 
light-cone momentum fractions and intrinsic momenta of all TMD distribution and fragmentation functions. On the 
contrary, all azimuthal dependences are explicitly shown. In particular, terms are collected according to the azimuthal 
dependence in the fragmentation process which directly enter the azimuthal asymmetries we want to study. Therefore, 
we get for the qq — > qq channel: 

[S^J + E^+vr)]™ 

-{/ qM A/s[I^i I 2 + I^2 I 2 + |m 3 | 2 ] 

-2A JV / a t M A JV / b t/ B M 2 M 3 °cos(^ 2 - <p 3 ) Jd*/, 
+ { A Ar / QTM / b/B M 1 °M 2 cos(^ - <p 2 ) 

-f a/A & N fbVBM°M° cos(^ - tps) } A N D n/q , cos $ . (21) 
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The symbol ~ is to recall that, as discussed above, this expression is valid only after integrating over the azimuthal 
angles of the initial intrinsic parton momenta, k± abl and is based on symmetry properties of the kernels under 
k±a,b —> — k±a,b- It contains less terms than the general expression for the kernels in the fcj^b-unintegrated, non 
planar partonic configuration. 

Analogously, for the numerator of the asymmetry, we find: 

~ { \A N f a/At f b/B cos0 a [ |M°| 2 + |M 2 °| 2 + |M 3 °| 2 ] 
-h a x A N f h y B cos(^ a - cp 2 + <^ 3 )2M 2 °M 3 
+ ^h^A N f byB cos(0 a + V2 - ^ 3 )2M 2 °Af 3 ° } smcj>s A D w/q 

+ { - ^A N f a/A ,A N f b y B cos0 a cos(^ 1 - <p 3 )M?M§ 

+Kfb/B COs(0 a +ipi~ <p 2 )M°M° 

-^h^f b/B cos(^ a + ^ 2 )M 1 °M 2 ° } sin^A^A^t cos^f 
+ { - \a n f a/A ,A N f b y B sin0 Q sm(Vi ~ ^)M°M 3 
-hlf b /B cos(0 a + <pi- y 2 )M 1 M 2 

-^h^h/B cos(0 a - ^ + ¥' 2 )M 1 °M 2 ° } coa<f>s A A N D K/9 i sin^f . (22) 

Let us discuss the physical content of these results. Eq. (|2"T|) gives the contribution of the — > channel to (twice) 
the unpolarized cross section. It contains two terms azimuthally symmetric in the fragmentation process: the first is 
the usual term already present in the collinear factorization scheme, the second one is the possible contribution due 
to the Boer-Mulders effect coming from both initial partons. 

The last two terms in Eq. (|2ip might potentially give rise to an azimuthal asymmetry in the jet — > tt + X unpolarized 
process: they are related to the combined action of the Boer-Mulders function (either for parton quark a or b 
separately) and of the Collins fragmentation function. Notice that only the first contribution in Eq. (|2"Tj) survives in 
a purely collinear scheme, or even in a scheme, like that adopted in Ref. 0, IH[, where intrinsic motion is kept 
into account only in the fragmentation process. 

Eq. (|22p is related to transverse spin asymmetries for pion production inside a jet. Again, it contains terms related 
to the unpolarized pion fragmentation function D 7T / q (z, k± 7T ), which are symmetric with respect to (j>^ , and terms 
proportional to the Collins fragmentation function which are responsible for the azimuthal asymmetries in the jet 
fragmentation process. 

The first group of terms, proportional to D^/q, are related to a single spin asymmetry (only hadron A is polarized). 
The only contribution allowed by rotational invariance and parity conservation comes from hadron A being polarized 
transversely to the jet production plane, which explains the sin<^s A term. 

Notice once more that the appearance of only the physically allowed contributions is not trivial in our expression, 
which, at this stage, is still unintegrated over fcj_ a ,6- However, as discussed above, taking into account symmetry 
properties under k± a-b — > — k± a ^ b amounts to select from the beginning only physically allowed contributions from 
the wealth of partonic terms present in the general non planar configuration for the partonic process. 

In the case of Eq. (|22[) the terms proportional to D^/ q come from the Sivers effect (first term) and from combinations 
of the transversity and the Boer-Mulders functions for partons a and b respectively. 

Let us now consider the terms in Eq. (|2"2")l related to the Collins fragmentation function, A N D 7r / q t. These refer 
effectively to a double spin asymmetry, since both hadron A and the final quark c (generating the observed jet) are 
transversely polarized. Their physical content is also easy to understand. The first three terms, proportional to 
A N D„i q t sm<f>s A , correspond to a double transverse spin asymmetry, where both hadron A and the final parton c 

are transversely polarized w.r.t. the jet production plane, along the Y cm axis. In this case, the Collins effect in the 
fragmentation process survives only for the component of the pion transverse momentum (w.r.t. the jet) orthogonal to 
the parton c polarization vector, that is the component laying in the production plane, which explains the associated 
cos0^ factor (see Fig. [T|). 

Analogously, the three terms proportional to A N D^^i cos tfis A correspond again to a double transverse spin asym- 
metry, this time for hadron A and the final parton c transversely polarized w.r.t. their own direction of motion but 
in the production plane (i.e. along the x axis of the respective helicity frames). Again, only the component of the 
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pion transverse momentum orthogonal to the parton c polarization vector contributes to the Collins asymmetry in 
the fragmentation process, which is this time guaranteed by the sin (j>^ factor. 

Although Eq. (f22j) makes the physical content of the asymmetry more evident, the following equivalent expression, 
where the terms proportional to the Collins FF are collected differently, makes the possible azimuthal asymmetries 
in the jet fragmentation process more readable: 



{ ±A N f a/A tf b/B cos0 a [ |M°| 2 + |M°| 



\M" 



~hlA N f b y B cos(0 a - tp 2 + if 3 )2M°M° 
+^h^A N f bf/B cos(0 a + ^ 2 - <p 3 )2M°M$ } sm<fi SA D„ 
{ [h a J b/B cos(0 a + ift - <p 2 )M°M° 



li 



--A N f a/A ,A N f b y B cos^ + ( p 1 - ^ 3 )M 1 °M 3 J sin(0 s , - <fg) 



r P 
* K ±a 



2M 



zhirh/B cos(^„ -tpi + <p 2 )M°M° 



A iV /a/At A N f b y B COs(0 a - i Pl + <p 3 )M»M» 



A N D 



■n/qf 



(23) 



Therefore, apart from the term proportional to D w / q , two azimuthal asymmetries in the distribution of leading 
pions inside the jet are possible, proportional respectively to sin(0g A =F <t>^)- 

Neglecting intrinsic motion of the initial partons, fc± ,& — > 0, Eqs. (|21[) . ([23| simplify considerably: 



[S(< 
[=( 



i s J + S(0 Sa + tt) ]™ f a/A (x a )f b/B (x b ) [ |M°| 2 + |M°| 2 + |M 3 T ]Ar/,(*, k ±7T ) , 
0S A ) - £(<fe A +7r)] w_ ' M ^(x a )/ b/s (x b )Af 1 °M 2 o A Ar J D 7r/9 t(z,fc±.)sin(0 SA - 0?) , 



(24) 
(25) 



in agreement with the results of Ref. |4£| . Notice however that our angle <p^ is the azimuthal angle of k ± w measured 
in the jet (parton c) helicity frame. Therefore, it does not coincide with the angle (f>h utilized in Ref. [40|, which is 
the azimuthal angle of k± v measured in the hadronic cm. frame [we call this angle 4>k in this paper, see Eq. ([7]) and 
Fig. H]. Only for forward jet production (cos#j — ¥ 1) do these angles coincide with good approximation. 
In the case of the gg — > gg partonic channel, the analogues of Eqs. (|2"Tj) . (|23|) are: 



[£(^ A ) + £(<fe, +tt) 

~ {fa/Afb/BWM^ + Wi' 



gg^gg 

0i2 , 



*/g 



~2A N f Ti y A A N f rib/B M°M$cos(v2 - ¥> 3 ) }Ar/ 
+ { A N f T « /A f b/B M°M° cosfai - If*) 

+f a/A A N f T y B M°M° cos(^ - cp 3 ) }A N D 7l/rf cos20f 

[E(0 Sa )-E(0 Sa +.)] 9 ™ 

~ { ~A N f a/A ,f b/B cos <t> a [ |M°| 2 + |M 2 °| 2 + |M 3 °| 2 ] 



(26) 



-lmF a +ZA"f Tib/B cos(0 Q - ^ 2 + ^ 3 )2M 2 U M 3 U 
-Im^+t A N f n/B cos(0 a + ip 2 - <^ 3 )2M 2 °Af 3 } sm</> 



S A Ar/<? 



1 



+-A v / a/At A^/tj/b cos(0 a - Vl + <p s )M?M$ sin(0 SA - 2$) 



lmF a - + tf b/B COs(0 a - if! + ^ 2 )M 1 °M 2 

1 



A N f a/Af A N f TNB cos(0 a + ^ - ^M^Mg ] sin(0 SA + 20f) Ja^D. 



t/7? 



(27) 
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The structure is the same as for the quark case, but this time all distributions related to transversely polarized 
quark partons and the Collins fragmentation functions are replaced by analogous functions for linearly polarized 
gluons, see Ref. [38| for more details. Notice that for linearly polarized gluons inside the polarized hadron A, to avoid 

confusion with notation we prefer to keep the definitions in terms of the functions F. °\? . It is also important to 

A A' A A 

notice that this time the possible azimuthal asymmetries in the distribution of leading pions inside the (gluon) jet are 
proportional to cos2</>^ and sin(4s A T 2</>^), respectively for the unpolarized and single polarized case. Therefore, 
by measuring these asymmetries one should in principle be able to select contributions coming from cither quark or 
gluon jet fragmentation. 

It is easy to see that in the case of collinear initial partons Eq. (|26[) reduces again to the usual collinear contribution 
to the unpolarized cross section, while Eq. (|27p vanishes. Therefore, the measurement of such types of asymmetries 
would be a clear indication that effects originating from intrinsic parton motion in the initial colliding hadrons are at 
work. From the phenomenological point of view, this would be a crucial test for the TMD approach, independently 
of the open issues concerning factorization and universality of the TMD distribution functions mentioned in the 
introduction. 

Expressions similar to those shown above for the qq — > qq and gg — > gg channels hold also for all partonic contri- 
butions involved, with the appropriate combinations of quark and gluon distribution and fragmentation functions. In 
general, less terms are present both in the denominator and the numerator of the asymmetry. Moreover, as a general 
rule quark (gluon) distribution and fragmentation functions off-diagonal in the parton helicity indexes, therefore as- 
sociated to transversely (linearly) polarized quarks(gluons), appear only in couple. This limits the number of allowed 
terms. 

According to these results, the single transverse polarized cross section for the process A(Sa)B — > jet + tt + X will 
have the following general structure: 

2da(4s A ,4%) ~ d a + dAa sin 4 Sa + doi cos <j% + dAof sin(0 s ^ ~4*) + dAcr+ sin(0 SA + $ ) 

+ dcr 2 cos 2$ + dAer^ sm(<f>s A -24>") + dAer+ sin(<f)s A + 24") . (28) 

Equivalently, the numerator and denominator of the asymmetry will have the following expression: 

da(4s A , 4") - da{<t> SA + tt, <j%) 

~ dAcr sin <f> SA + dAcr~ sin(4 SA -(/>") + dAcr+ sm(4 SA + 

+ dAa 2 " sin(0 SA -24") + dAa+ sin(0 SA +2$), (29) 



M4s A ,<$) + da(4s A +7T,4") = 2da un P(0f ) ~ da + da x cos $ + da 2 cos 2<fi . (30) 

In terms of the polarized cross section, Eq. (|28p , we can define average values of appropriate circular functions of 
(j>s A and (j>„ , in order to single out the different contributions of interest: 

/w( , , H ^ f . x [Ws^l W(4s A ,4")da(4s A ,4") 

( w(4 SA ,4,))(p,^) = Sd4> s A <WM*s A ,ti<) ■ (31) 

Alternatively, for the single spin asymmetry we can, in close analogy with the case of semi-inclusive deeply inelastic 
scattering, define appropriate azimuthal moments, 

A^ 3 ^\p h z,k^) = 2(W(4s A ,^))(p- 3 ,z,k^) 

= J d4s A d^ W(4s A , 4? ) [da(4s A , <$) ~ da(4s A + tt, )] 

fd4 SA d4%[da(4s A ,4%) + da(4s A +TT,4%)} ' 1 ] 

where W(4s A j 4") i s again some appropriate circular function of 4s A and 4" ■ I R practice, it will be any of the 
circular functions appearing e.g. in Eqs. (|2"3")l , (|2T|) for specific partonic channels, and for polarized cross sections in 
general in Eq. (f2T)f so that the coefficient related to the corresponding azimuthal moment is singled out. 



III. PHENOMENOLOGY 



In this section we will present and discuss some phenomenological implications of our approach for the unpolarized 
and single-transverse polarized case in kinematical configurations accessible at RHIC by the STAR and PHENIX 
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experiments. We will consider both central (?jj = 0) and forward (rj s = 3.3) (pseudo)rapidity configurations and 
different cm. energies, y/s = 62 A, 200, 500 GeV, aiming at a check of the potentiality of the approach in disentangling 
among different quark and gluon originating effects. We will also consider two very different situations as far as concern 
the TMD distribution and fragmentation functions involved. 

We will first consider, for ir + production only, a scenario in which the effects of all TMD functions are over- 
maximized. By this we mean that all TMD functions are maximized in size by imposing natural positivity bounds 
(and the Soffer bound for transversity [461 \4l\); moreover, the relative signs of all active partonic contributions are 
chosen so that they sum up additively. This very extreme scenario of course might imply the violation of other, more 
stringent, bounds and sum rules; examples are the Burkardt sume rule for the Sivers distribution (48|, and the Schafer- 
Teryaev sum rule for the Collins function 49] . On the other hand, it has the advantage of setting an upper bound on 
the absolute value of any of the effects playing a potential role in the azimuthal asymmetries. Therefore, all effects 
that are negligible or even marginal in this scenario may be directly discarded in subsequent refined phcnomcnological 
analyses. 

As a second step in our study we will consider, for both neutral and charged pions, only the surviving effects, 
involving TMD functions for which parameterizations are available from independent fits to other spin and azimuthal 
asymmetries data in SIDIS, DY, and e + e~ processes. Although in our approach factorization and universality are 
not guaranteed for the process under consideration, we still believe that at the present stage this analysis can be 
of phcnomenological relevance. It can certainly help in pointing out inconsistencies among fits based on different 
processes, that could be a signal of universality-breaking effects. On the contrary, good consistency among fits to 
different observables from SIDIS, e + e~ and hadronic collisions data, while not proving factorization, might signal 
the smallness of possible universality breaking terms and the usefulness of the factorization hypothesis in the present 
phenomenological analyses. 

In this paper, for numerical calculations all TMD distribution and fragmentation functions will be taken in the 
simplified form where the functional dependences on the parton light-cone momentum fraction and on transverse 
motion are completely factorized, assuming a Gaussian-like flavour-independent shape for the transverse momentum 
component. Preliminary lattice QCD calculations seem to support the validity of this assumption, see e.g. Ref. [Hoj|. 
Notice however that kinematical cuts introduced to prevent, as usual in the parton model, that the parton longitudinal 
momentum (energy) be opposite to (larger than) that of the parent hadron, effectively lead to a correlation between 
the light-cone momentum fraction and the transverse momentum, particularly at very small and very large (— > 1) 
light-cone momentum fractions (see e.g. appendix A of Ref. (4l|V 

For the generic parton a, the unpolarized and any of the polarized functions (that is, the Sivers and Boer-Mulders 
distributions and the Collins FF, and the analogous ones for gluons) will therefore assume respectively the form 
T™v(u,p) = / Q unp (u)3 (0) (p) and AT a (u,p) = Af a (u)g^(p), with a = q,q,g, i = 0, 1, 2, 3, u = x or z and p = k± 
or kj_n for distribution/fragmentation functions respectively. Our parameterizations are required to respect angular 
momentum conservation in the forward direction, Eqs. (|lip . (1171) ; therefore we define 

ff W (p)=(^)^ W (p)- (33) 
In particular, g^°\p) = g(p) is a simple Gaussian normalized to unity: 

3<°>(p) = g(p) = exp[ -p 2 /{pD] , (34) 

while 

9 {l) (p) = (^y hMHp)^K t (J^y exp[-p 2 /<^>], 1 = 1,2,3. (35) 

All the polarized TMD functions are required to fulfill natural positivity bounds (for transversity, the Soffer bound) 
with respect to the correponding unpolarized functions, coming from their general definition as 

F(S) + F(-S) nJ^P ' 1 ' 

where S is here the spin of the polarized quark or hadron involved, and n = 1,2, depending wether the polarized 
particle is respectively the final/initial one in the soft process considered (analogous relations hold for gluons). The 
positivity bound therefore reads: 



(37) 
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As a matter of fact, to simplify relations we will consider a more conservative and stringent bound on the two factored 
components, 

nf a p (m) 9 (0, (P) 
The first condition is usually fulfilled by defining 

Af a (u) = nN a (u)fr P (u) = nN a u a °(l - uf« { a + P - a> , with |JV„| < 1 . (39) 

Oi a a Pa 

Concerning the transverse momentum dependent component, g^\p), the positivity bound, Eq. (|38[) . can only be 
fulfilled if, in Eq. (|35|) , (pf) < (pq). We then fix the factors K t in Eq. (|35)l by saturating the bound at the maximum 
value of Finally, once a choice has been performed for (pg ), the (pf ) are fixed by maximizing the corresponding 
(i + l)-th p-moments. It is then an easy exercise to verify that this gives the conditions: 

<*?> = §<*?>, (A) = \<P&), (pI) = 1(pD- (40) 

From the above equations it is clear that the over-maximized scenario for 7r + production we are going to present 
is obtained by taking, for all polarized TMD functions, the coefficients Af a (u) = 1. For the transverse momentum 
component, the procedure delineated above guarantees the correct power-like behaviour in the forward direction, 
while maximizing the moments of the involved functions. 

One also needs LO parameterizations for the usual unpolarized collinear distribution and fragmentation functions. 
Concerning the parton distribution functions, we will adopt the unpolarized set GRV98 [5l[ and (for the Soffer bound) 
the corresponding longitudinally polarized set GRSV2000 [52]. Since the range of the jet transverse momentum (the 
hard scale) covered is significant, we will take into account proper evolution with scale. Concerning transversity, in 
the maximized scenario we will fix it at the initial scale by saturating the Soffer bound and then letting it evolve. 
On the other hand, the transverse momentum component of all TMD functions is kept fixed with no evolution with 
scale. Notice that at this stage evolution properties of the full TMD functions are not known. 

As for fragmentation functions, we will adopt two well-known LO sets among those available in the literature, the 
set by Kretzer (K) [53[ and the one by De Florian, Sassot and Stratmann (DSS) [54| . Our choice is dictated by the 
subsequent use of the two available parametrization sets for the Sivers and Collins functions in our scheme, that have 
been derived in the past years by adopting these sets of FFs. Let us notice that, as the authors of Ref. [5J] suggest, 
the LO DSS fragmentation function set has to be handled with some care. In fact, aiming at reproducing for the 
first time unpolarized cross sections for inclusive hadronic collisions, at LO accuracy a very huge gluon component 
(as compared to other FF sets) is required. This could be an artifact of the LO set, which in fact is sizably reduced 
in the NLO parameterizations. However, we are at present forced to work at leading order in our TMD approach. 
The comparison among the two sets adopted will allow to stress the possible effects of the large gluon component in 
the LO DSS set. 

Concerning the parameterizations of the transversity and Sivers distributions, and of the Collins functions, we will 
consider two sets resulting from the fits to available data on azimuthal asymmetries in polarized SIDIS from HERMES 
and COMPASS experiments, and on hadron-pair production in e + e~ collisions from Belle: 

Set 1 (SIDIS 1) includes the u, d quark Sivers functions of Ref. [55| . the u, d quark transversity distributions and 
the favoured and unfavoured Collins FFs of Ref. [56[ . Data include preliminary HERMES data for charged pions [13] 
and COMPASS data on charged hadrons with a deuteron target [58| on the SIDIS Sivers asymmetry; HERMES data 
for charged pions [IJ H^l an d COMPASS data for charged hadrons with a deuteron target [6(| on the SIDIS Collins 
asymmetry; early Belle data on azimuthal asymmetries for hadron-pair production in e + e~ collisions [14]; no SIDIS 
data on kaons were used and the Kretzer set [53[ for pion FFs was used. 

Set 2 (SIDIS 2) includes the new u, d, and sea-quark Sivers functions of Ref. [6lj and an updated set of the it, d quark 
transversity distributions and of the favoured and unfavoured Collins FFs of Ref. [62j. Corresponding data include: 
preliminary HERMES data on pions and charged kaons [63[ and preliminary COMPASS data on charged pions and 
kaons with a deuteron target [64| for the SIDIS Sivers asymmetry; preliminary HERMES data for pions [63| and 
COMPASS data for charged pions with a deuteron target [l2j on the SIDIS Collins asymmetry; recent Belle data on 
azimuthal asymmetries for hadron-pair production in e + e~ collisions (l5| : the DSS set |54j for pion and kaon FFs was 
adopted. 

Notice that the almost unknown gluon Sivers function was tentatively taken positive and saturated to an updated 
version of the bound obtained in Ref. [65[ by considering PHENIX data for the tt° transverse SSA at mid-rapidity 
production in polarized pp collisions at RHIC (66j j . 
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It is also important to stress here that polarized SIDIS data on azimuthal asymmetries from HERMES and COM- 
PASS experiments cover a relatively limited range of Bjorken x, x^ < 0.3. Therefore, the statistical uncertainty 
of the parameterizations available for the transversity and Sivers distributions is huge at large x values, where one 
extrapolates their behaviour. As we will see, this reflects in the very different behaviour of the Sivers and Collins 
asymmetries when estimated adopting sets SIDIS 1 and SIDIS 2 respectively. 

Finally, regarding the quark Boer-Mulders distribution function much less is known and available parameterizations 
have large uncertainties. In our calculations we have adopted the recent parametrization by Barone, Melis and 
Prokudin (BMP) 67], which makes use of our Set SIDIS 2 for the transversity and Sivers distributions and for the 
Collins function. 

We have considered the following kinematical configurations for the PHENIX and STAR experiments at RHIC: 

1) = 62.4 GeV, rjj = 0, 1 < p jT < 14 GeV; 

2) At yfs = 62.4 GeV, the forward rapidity configuration covers a very limited range of p^x values which probably 
prevents the unambiguous definition of jets and leading particles, therefore we will not consider it in the sequel; 

3) V« = 200 GeV, rjj = 0, 2 < p ]T < 15 GeV; 

4) y/s = 200 GeV, = 3.3, 2 < p ]T < 6.5 GeV [0.27 < x F < 0.88]; 

5) = 500 GeV, rjj = 0, 2 < p ]T < 15 GeV; 

6) y/s = 500 GeV, t/j = 3.3, 2 <p iT < 15 GeV [0.11 < x F < 0.81]. 

For completeness, in the forward rapidity case we have also shown the range of x F covered, where x F is the usual 
Feynman variable for the jet, x F — 2p j i J /y/s. 

In all cases considered, since we are interested to azimuthal asymmetries for leading particles inside the jet, we will 
present results obtained integrating the light-cone momentum fraction of the observed hadron, z, in the range z > 0.3. 

A. Azimuthal asymmetries for the unpolarized cross section 

In this section we will discuss results for the azimuthal (cos0^), (cos2</>^) asymmetries [see Eq. pip ] in the 
unpolarized cross section for the process pp — > jet + tt + X. 

As it is easy to verify by looking e.g. at Eqs. (f2"T|) . (|2^|) for the qq — » qq and gg — > gg partonic contributions 
(analogous results, where allowed, hold for all other channels), in the unpolarized case: 

1) The symmetric part gets contributions by the usual unpolarized term, already present in the collinear approach, and 
by an additional term involving a Boer-Mulders £g> Boer-Mulders convolution for the initial quarks (or the analogous 
terms involving linearly polarized gluons); however, we have explicitely checked that even in the maximized scenario 
this last contribution is always negligible in all the kinematical configurations considered; therefore, we will not discuss 
it anymore in the sequel; 

2) The cos <f>^ asymmetry is generated by the quark Boer-Mulders £g) Collins convolution term, involving a transversely 
polarized quark and an unpolarized hadron both in the initial state and in the fragmentation process. In the central 
rapidity region {rjj = 0) the maximized value of this asymmetry is of the order 1-3%, depending on the fragmentation 
function set adopted and on the cm. energy considered, being almost negligible at y/s = 500 GeV. In the forward 
rapidity region, ?yj = 3.3, the maximized cos cj)^ asymmetry can be much larger both at y/s = 200 and 500 GeV. As 
an example, in Fig. [2] we show the maximized cos cf>^ asymmetry (solid red lines) for ir + production at cm. energy 
y/s = 200 GeV in the central (left panel) and forward (right panel) rapidity region as a function of p\t-, from pjt = 2 
GeV up to the maximum allowed value, adopting the Kretzer FF set. Slightly lower values are obtained using the 
DSSset. 

3) The cos 2c/)^ asymmetry is related to the term involving linearly polarized gluons and unpolarized hadrons both in 
the initial state and in the fragmentation process, that is the convolution of a Boer-Mulders-like gluon distribution with 
a Collins-like gluon FF. Even the maximized contribution is practically negligible in the kinematical configurations 
considered. As an example, again in Fig. [21 we show the maximized cos 2(f>^ asymmetry (dashed green lines) for ir + 
production at t/s = 200 GeV cm. energy in the central (left panel) and forward (right panel) rapidity region as a 
function of PjT, adopting the Kretzer FF set. Similar results are obtained using the DSS set. 

Concerning results with available parameterizations, for the quark-originated cos0^ asymmetry we have verified 
that the asymmetries obtained with the parameterizations adopted here, our set SIDIS 2 and the BMP set for the 
Boer-Mulders function, are negligible in all kinematical configurations considered. No parameterizations are presently 
available for the analogous gluon contributions leading to the cos 2</>^ asymmetry. 
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FIG. 2: (color online). Maximized quark-originated (cos (j>^ ) and gluon-originated (cos 2(f>^ ) asymmetries (solid-red and dashed- 
green lines respectively) for the unpolarized pp — > jet + 7r + + X process, at yfs = 200 GeV cm. energy in the central (left panel) 
and forward (right panel) rapidity region as a function of p\t, from pjr — 2 GeV up to the maximum allowed value, adopting 
the Kretzer FF set. Slightly lower (similar) values are obtained for quark(gluon) asymmetries when using the DSS set. 



B. Azimuthal asymmetries for An(p^P — > jet + it + X) 



Let us now discuss our numerical results for the Sivers (A N A ) asymmetry and the quark [A N A " ] and 

gluon [^4™ < -^ s a t2< ^t j Collins (-like) asymmetries, see Eq. ([32]). Our estimates are qualitatively similar at the three 
different cm. energies considered, with some differences in the size of the asymmetries and in the relative weight of 
the quark and gluon contributions where both play a role. Therefore, we will concentrate on the results obtained at 
y/s = 200 GeV. 



1. The Sivers asymmetry 

In this case, both quark and gluon contributions can be present, and they cannot be disentangled. However, some 
kinematical configurations can be dominated by quark or gluons terms and a sizable asymmetry in these regions might 
be an unambiguous indication for a Sivers asymmetry generated by the dominant partonic contribution. 

In Fig. [3] we show the total observable Sivers asymmetry (solid red line), and the corresponding quark and gluon 
contributions (dashed green and dotted blue lines respectively) for tt + production, in the maximized scenario and 
adopting the Kretzer fragmentation function set, at yfs = 200 GeV and as a function of pjr in the central (left 
panel) and forward (right panel) rapidity regions. The maximized potential Sivers asymmetry can be very large in 
both cases. In the central rapidity region, the asymmetry is dominated by the gluon contribution at the lowest pj t 
range while gets comparable quark and gluon contributions in the large pjT range. A large Sivers asymmetry around 
Pj t = 4 -i- 6 GeV could then be a clear indication for a sizable gluon contribution. However, one must not forget that, 
as mentioned above, recent PHENIX results for A^(p^p — > tt° + X) in the central rapidity region put much more 
stringent bounds on the gluon Sivers distribution than the simple positivity bound adopted in the maximized scenario. 
We have checked that even adopting this more stringent bound, a potential gluon-generated Sivers asymmetry of the 
order of 2% might survive in this region, being possibly measurable. In the forward rapidity region, on the contrary, 
the quark and gluon contributions are comparable at low pjT values, while the maximized asymmetry is dominated 
by the quark contribution for p } t ft 4 GeV. Therefore, a large Sivers asymmetry in this kinematical range could be 
ascribed unambiguosly to the quark Sivers effect. Qualitatively similar results are obtained at the other cm. energies 
considered or adopting the DSS set of fragmentation functions, with some changes in the total size and in the relative 
weights of the quark and gluon contributions. 

In Fig. 2] we show, for both neutral and charged pions, the quark and gluon contributions to the Sivers asymmetry, 
obtained adopting respectively the parametrization sets SIDIS 1 (quark contribution: solid red line; gluon contribution: 
dashed green line) and SIDIS 2 (quark contribution: dotted blue line; gluon contribution: dot-dashed cyan line), and 
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FIG. 3: (color online). Maximized total (solid red line), quark-originated (dashed green line) and gluon-originated (dotted blue 
line) Sivers asymmetry for the p^p — > jet + 7r + +X process, at yfs = 200 GeV cm. energy in the central (left panel) and forward 
(right panel) rapidity region as a function of PjT, from pjT = 2 GeV up to the maximum allowed value, adopting the Kretzer 
FF set. Similar results, with some differences in the total size and in the relative weight of the quark and gluon contributions 
are obtained adopting the DSS set of fragmentation functions and considering different cm. energies. 



the updated version of the bound on the gluon Sivers asymmetry derived in Ref. [65j, at y/s = 200 GeV and in the 
forward rapidity region, as a function of PjT- The dotted black vertical line delimits the region beyond which the 
SIDIS parameterizations for the quark Sivers distribution are extrapolated outside the xb region covered by SIDIS 
data and are therefore plagued by large uncertainties. This reflects on the fact that below this limit the two sets give 
comparable results, while above it they differ remarkably. In particular, at the largest reachable pjT values the SIDIS 1 
set gives a Sivers asymmetry of the order 2 4%, while the SIDIS 2 set leads to a negligible asymmetry. Therefore, 
a measurement of this asymmetry might help in clarifying the behaviour of the quark Sivers distribution in the large 
x region, which is not covered by present SIDIS data from HERMES and COMPASS experiments. Future planned 
measurements at the Jefferson Lab (Jlab) 12 GeV Upgrade will also be very useful in this respect (see e.g. Ref. [68| ) . 
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FIG. 4: (color online). The estimated quark and gluon contributions to the Sivers asymmetry for the p^p — > jet + 7r + X process, 
obtained adopting respectively the parametrization sets SIDIS 1 (quark contribution: solid red line; gluon contribution: dashed 
green line) and SIDIS 2 (quark contribution: dotted blue line; gluon contribution: dot-dashed cyan line), at ^/s = 200 GeV 
cm. energy in the forward rapidity region and as a function of PjT, from pjT = 2 GeV up to the maximum allowed value. 
The dotted black vertical line delimits the region beyond which the SIDIS parameterizations for the quark Sivers function are 
presently plagued by large uncertainties. Similar results are obtained when considering different cm. energies. 
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FIG. 5: (color online). Maximized quark (solid red line) and gluon (dashed green line) Collins(-like) asymmetries for the 
p^p jet + 7r + + X process, at ^/s — 200 GeV cm. energy in the central (left panel) and forward (right panel) rapidity region 
as a function of pj t, from pj t — 2 GeV up to the maximum allowed value, adopting the Kretzer FF set. Notice the difference in 
the scale between the two panels. Similar results, with some differences in the total size and in the relative weight of the quark 
and gluon contributions are obtained adopting the DSS set of fragmentation functions and considering different cm. energies. 



2. The Collins A N ^ s a^^^ ) an ^ Collins-like A N ^ Sa ^ ' asymmetries 

Let us first briefly discuss the quark generated asymmetry ^4™^ s ^ + ^ , ' r _ comes from two distinct contributions, 
see e.g. Eq. ([23]) for the qq — > qq channel: one involving the convolution between the term of the TMD transversity 
distribution suppressed in the collinear configuration (oc fcj^iT) and the Collins function; another term involving 
the convolution of the Sivers and Boer-Mulders distributions for the initial quarks with the Collins function for the 

final quark [an analogous term appears also in the A s ^^ Sa ' asymmetry, see Eq. (|2"3"1) ]. We have explicitly checked 
that for the process under study and the kinematical configurations considered both these contributions are always 
negligible already in the maximized scenario. Therefore we will not consider the s'm(cj)g A + <p^ ) asymmetry in the 

sequel. A similar situation holds also for the gluon generated j^™^ s a +2< ^^ ' asymmetry, where two contributions 
analogous to the quark ones discussed above but for linearly polarized gluons are involved. 

In Fig.[5]we present the quark A s ^^ Sa ^ ' Collins asymmetry (solid red lines) and the gluon A s ^^ Sa 2{ ^" ^ Collins- 
like asymmetry (dashed green lines) in the maximized scenario for the p^p — > jet + ir + + X process, at y/s = 200 GeV 
cm. energy in the central (left panel) and forward (right panel) rapidity region as a function of PjT, from pjT = 2 
GeV up to the maximum allowed value, adopting the Kretzer FF set. In the central rapidity region the quark Collins 
asymmetry is very small at the lowest p 3 t values, then increases almost linearly reaching about 8% at the upper 
range. At y^s = 500(62.4) GeV cm. energy the behaviour is similar and the largest reached size, at large p 3 T values, 
is about half (twice) respectively. Results with the DSS fragmentation function set are slightly lower in size. Instead, 
in the forward rapidity region the asymmetry is (potentially) always large and increases almost linearly from about 
25% to about 70% going from the lowest to the largest pj t values. Results are very similar at y/s — 500 GeV and 
when adopting the DSS fragmentation function set. 

Concerning the gluon Collins-like asymmetry, both in the central and in the forward rapidity regions it is of the 
order of 5% at the lowest p } t values, then starts decreasing slowly and becomes negligible at large pj t values. Very 
similar results hold at different energies and when adopting the DSS set. 

Let us now consider, for both neutral and charged pions, numerical results for the quark Collins asymmetry obtained 
adopting the parameterizations SIDIS 1 and SIDIS 2 for the transversity distribution and the Collins fragmentation 
function (no parameterizations are available yet in the analogous gluon case) . It turns out that in the central rapidity 
region the estimated asymmetry is practically negligible in all cases considered (different cm. energies and FF sets). 
Only for the SIDIS 2 parametrization and at = 62.4 GeV the asymmetry for charged pions can reach about 2-^3% 
in size at large p- } T values. 

Concerning the forward rapidity region, in Fig. [6] we present, for both charged and neutral pions, some results at 
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FIG. 6: (color online). The estimated quark Collins asymmetry for the p^p — > jet + it + X process, obtained adopting the 
parameterizations SIDIS 1 (left panel) and SIDIS 2 (right panel) respectively, at i/s — 200 GeV cm. energy in the forward 
rapidity region and as a function of PjT, from p^r = 2 GeV up to the maximum allowed value. Notice the difference in the 
scale between the two panels. The dotted black vertical line delimits the region beyond which the SIDIS parameterizations for 
the quark transversity distribution are presently plagued by large uncertainties. Similar results are obtained when considering 
different cm. energies. 

= 200 GeV cm. energy as a function of pjT, adopting the SIDIS 1 (left panel) and the SIDIS 2 (right panel) 
parameterizations. The Collins asymmetry for neutral pions (dashed green lines) results to be practically negligible. 
This can be easily understood since in the available parameterizations the favoured (e.g. u — > 7r + ) and unfavoured 
(e.g. d — > 7r + ) Collins fragmentation functions are comparable in size and opposite in sign. Due to isospin symmetry 
the 7T° FFs are half the sum of those for charged pions, therefore the tt° Collins FF is always very small. Moreover, 
the u, d quark transversity distributions are also opposite in sign, leading to additional cancellations among quark 
contributions. For charged pions, similarly to the case of the Sivers asymmetry, the two parameterizations give 
comparable results (notice the different scale adopted in the two panels) in the pj t region where the transversity 
distribution is reasonably constrained by SIDIS data (see the dotted black vertical line), while they lead to completely 
different estimates in the large p- 3 t region where the parameterizations are basically unconstrained by SIDIS data. In 
particular, in this region the SIDIS 1 set gives almost negligible results, while the SIDIS 2 set leads to an asymmetry 
of about 8% (tt + ) and 15% (tt _ ) in size, which should be hopefully measurable. A measurement of this asymmetry 
would be then very important and helpful in clarifying the large x behaviour of the quark transversity distribution. 
A qualitatively similar situation is obtained at ^fs = 500 GeV. 

3. Transverse single spin asymmetry for inclusive jet production 

For completeness we have extended our analysis to the transverse single spin asymmetry for inclusive jet production 
in polarized pp collisions, A^{p^p — > jet + X). In principle, this case can be obtained by the jet + pion production 
process by integrating over the full pion phase space. Of course, in this case the unobserved fragmentation process 
in the final state plays no role in the azimuthal asymmetries, which can only be originated by mechanisms, like the 
Sivers effect, acting in the initial state. Moreover, we have verified that for the kinematical configurations considered 
in this paper all contributions but the Sivers effect play a negligible role already in the maximized scenario. Therefore, 
in what follows, we limit our discussion to the Sivers asymmetry. As already mentioned, in this case quark and gluon 
contributions cannot be disentangled since they add up leading to a sin (j>s A asymmetry. 

Let us first discuss the maximized scenario. In the central rapidity region, the maximized gluon contribution is 
of the order 20% at the lowest pj t values, decreasing fast to about 3% at large p 3 - t for all cm. energies considered. 
The maximized quark contribution is of the order 1-3% in the full p 3 t range, slowly decreasing with the increase of 
the cm. energy. The total potential effect is therefore sizable only at small p 3 T values due to the gluon component. 
The situation is different in the forward rapidity region. Here both quark and gluon maximized contributions can be 
very sizable, showing as expected an opposite, respectively increasing and decreasing, behaviour vs. p 3 T- The total 
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FIG. 7: (color online). The estimated quark and gluon Sivers contributions to the transverse single spin asymmetry for the 
p^p — > jet + X process, at y/s = 200 GeV cm. energy in the central (left panel) and forward (right panel) rapidity region as 
a function of PjT, from pjT = 2 GeV up to the maximum allowed value, obtained adopting the parametrization sets SIDIS 1 
(quark contribution: solid red line; gluon contribution: dashed green line) and SIDIS 2 (quark contribution: dotted blue line; 
gluon contribution: dot-dashed cyan line). The dotted black vertical line in the right panel delimits the region beyond which 
the SIDIS parameterizations for the quark transversity distribution are presently plagued by large uncertainties. Similar results, 
with some differences in the total size and in the relative weight of the quark and gluon contributions are obtained considering 
different cm. energies. 



maximized Sivers effect is therefore large in the full pj t range with little dependence on the cm. energy. 

Concerning numerical estimates obtained adopting the available parameterizations SIDIS 1 and SIDIS 2 for the 
quark Sivers function, and the updated bound on the gluon Sivers function, the situation is the following: 

1) In the central rapidity region, for both SIDIS 1,2 sets and all energies considered the quark contribution is practically 
negligible. Instead, the gluon contribution can be at most of the order 10 15% at the lowest pj t values but decreases 
quickly with the increasing of pjx- However, at least for y/s — 200 and 500 GeV, it can still be about 2-4% in the 
upper pj t range. The measurement of a comparable Sivers asymmetry in these kinematical configurations could then 
be a clear indication for a gluonic contribution to the Sivers effect. 

2) In the forward rapidity region the quark contribution is small and negative at pj t = 2 GeV for both sets adopted, 
while at large p } T values it is negligible for the SIDIS 2 set and positive and of the order 2-4% for the SIDIS 1 set. 
The gluon contribution can be sizable at very low pjT values but becomes negligible quickly as pjT increases. 

As an example, in Fig. [7] we show the estimated quark and gluon Sivers contributions to the transverse single 
spin asymmetry for inclusive jet production in the central (left panel) and forward (right panel) rapidity regions at 
y/s = 200 GeV, obtained adopting the SIDIS 1 and SIDIS 2 parameterizations for the quark Sivers function and the 
updated bound for the gluon Sivers function (assumed to be positive). 



IV. CONCLUSIONS 



In this paper we have presented a study of the azimuthal asymmetries measurable in the distribution of leading 
unpolarized or spinless hadrons (mainly pions) inside a large-j>r jet produced in unpolarized and single-transverse 
polarized proton proton collisions for kinematical configurations accessible at RHIC. To this end, we have adopted 
a generalized TMD parton model approach with inclusion of spin and intrinsic parton motion effects both in the 
distribution and in the fragmentation sectors. We have shown how a detailed phenomenological analysis of these 
effects can be very useful in shedding light on several aspects of azimuthal and transverse single spin asymmetries in 
(un)polarized hadronic collisions. It may also help in clarifying the role played by the quark(gluon) Sivers distribution 
and by the Collins (-like) fragmentation function in the sizable single spin asymmetries observed at RHIC for forward 
pion production. Available parameterizations for the TMD quark transversity and Sivers distribution functions, 
obtained by fitting SIDIS and e + e~ data, are presently largely unconstrained for light-cone momentum fractions 
x > 0.3, that is the region playing a fundamental role for forward pion production at RHIC. The transverse single spin 
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asymmetry for inclusive particle production is a complicate higher-twist effect involving several TMD mechanisms 
that cannot be easily disentangled as in the case of SIDIS and DY processes. On the contrary, the leading-twist 
azimuthal asymmetries discussed in this paper allow one, in close analogy with the SIDIS case, to discriminate among 
different effects by taking suitable moments of the asymmetries. Moreover, we have shown that in principle quark 
and gluon originating jets can be distinguished, at least in some kinematical regimes. Neglecting intrinsic motion in 
the distribution sector leaves at work only the Collins azimuthal asymmetry. As already shown by Yuan [40l |. the 
measurements proposed in this paper would allow one to determine unambiguously the role played by the Collins 
effect. Universality properties of the Collins function have been proved, so that this process can be complementary 
to the SIDIS and e + e _ measurements in order to constrain the quark Collins fragmentation function and, as an 
important by-product, the large x behaviour of the TMD transversity distribution function. 

For the quark and gluon Sivers function and for the Boer-Mulders function the situation is more complicate. Since 
factorization has not been proven yet for inclusive particle production in hadronic collisions, the use, based on univer- 
sality, of the parameterizations obtained by fitting SIDIS and e + e~data might be questionable. The phenomenological 
analysis proposed in this paper gives us the opportunity of testing the factorization and universality assumptions, and 
of gaining information on the size and sign of the TMD functions discussed. This can be very useful also for further 
developments of the TMD approach, since it is difficult a priori to assess at which values of the factorization scale the 
role and the size of possible factorization-breaking terms become relevant and non negligible. 

Let us finally stress again that the unambiguous measurement of any of the asymmetries, other than the Collins 
one, discussed in this paper would be a clear indication of the role played by intrinsic parton motion in the initial 
colliding hadrons for the spin asymmetry sector in polarized hadronic collisions. 
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